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Abstract 

An  efficient  and  powerful  computational  method  is  presented  to  synthesize  a  meta-material  to  specified  electro¬ 
magnetic  properties.  Using  the  periodicity  of  meta-materials,  the  Finite  Element  Methodology  (FEM)  is  developed  to 
estimate  the  reflection  and  transmission  through  the  meta-material  structure  for  a  normal  plane  wave  incidence.  For 
efficient  computations  of  the  reflection  and  transmission  over  a  wide  band  frequency  range  through  a  meta-material 
a  Finite  Difference  Time  Domain  (FDTD)  approach  is  also  developed.  Using  the  Nicholson-Ross  method  and  the 
Genetic  Algorithms,  a  robust  procedure  to  extract  electromagnetic  properties  of  meta-material  from  the  knowledge 
of  its  reflection  and  transmission  coefficients  is  described.  Few  numerical  examples  are  also  presented  to  validate  the 
present  approach. 


Introduction 

In  the  recent  years  there  has  been  a  great  deal  interest  in  studying  the  material  (called  meta-material  [1-3])  hav¬ 
ing  simultaneously  negative  permittivity  and  permeability.  These  meta-materials  do  not  exist  in  nature  and  are  artifi¬ 
cially  formed  by  embedding  properly  designed  metallic/non-metallic  inclusions  in  binding  mediums.  In  such  a 
)  A 

material,  the  wave  vector  k  and  the  Poynting  vector  P  are  anti-parallel  causing  reversal  of  some  basic  electromag¬ 
netic  wave  propagation  characteristics  such  as  Snell’s  law  and  Doppler  effects.  The  reversal  of  these  basic  propaga¬ 
tion  properties  has  lead  to  some  potential  applications  of  these  materials  such  as  high  resolution  lenses, 
sub  wavelength  resonators  and  enhancement  of  antenna  gains.  One  of  the  challenges  faced  by  the  researchers  in  this 
field  is  lack  of  a  robust  method  to  characterize  these  materials.  The  characterization  of  the  meta-materials  is  carried 
out  in  two  steps.  First,  the  reflection  and  transmission  properties  of  these  materials  are  estimated  /measured  using  a 
computational/experimental  method.  In  the  second  step,  a  proper  inversion  process  is  used  to  extract  effective  elec¬ 
tromagnetic  properties  from  the  knowledge  of  reflection  and  transmission  coefficients.  Experimental  procedures  are 
too  expensive  and  time  consuming  to  explore  all  design  space,  and  are  is  usually  not  preferred  in  early  stages  of 
design  cycle.  In  the  present  paper,  an  efficient  and  powerful  computational  method  is  presented  to  characterize/esti¬ 
mate  effective  electromagnetic  properties  of  meta- materials.  In  this  method,  first,  the  metallic  inclusions  of  interest 
are  modeled  using  a  commercially  available  Computer  Aided  Design  (CAD)  package.  Using  the  Finite  Difference 
Time  Domain  method  the  reflection  and  transmission  through  these  inclusions  are  then  estimated.  In  the  final  step 
using  an  inversion  procedure  based  on  two  approaches,  Nicolson-Ross  method  [4]  and  Genetic  Algorithm,  effective 
properties  of  the  selected  inclusions  are  extracted.  Few  examples  with  known  effective  properties  are  presented  to 
validate  the  present  method. 


Analysis  of  Meta-Materials 

In  this  section,  the  development  of  computational  methods  based  on  the  Finite  Element  Methodology  and  also 
Finite  Difference  Time  Domain  (FDTD)  method  are  presented.  Figure  1(a)  shows  a  metallic  inclusion  in  a  split  ring 
form  and  Figure  1(b)  shows  a  dielectric  medium  embedded  with  a  periodic  arrangement  of  split  ring  shown  in  Figure 
1(a).  It  is  assumed  that  the  structure  is  infinite  in  the  x-y  plane  and  has  finite  width  in  the  z-direction.  It  also  assumed 
that  the  structure  is  illuminated  by  a  plane  wave  with  normal  incidence  as  shown  in  Figure  1(b). 
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Figure  1(a)  Split  ring. 


J - j/ — Split  Rings 

Figure  1(b)  Dielctric  slab  embedded  with  split  rings. 


Finite  Element  Method:  To  estimate  the  reflection  and  transmission  coefficients  of  such  a  structure,  a 
unit  periodic  cell  as  shown  in  Figure  1(b)  is  considered.  Using  the  electric  and  magnetic  symmetries,  the 
unit  cell  can  be  terminated  by  electric  and  magnetic  conducting  boundaries  at  y  =  +a/2  and 
x  =  +b/2,  respectively.  The  electromagnetic  field  incident  at  the  input  port  of  the  unit  cell  can  be  writ¬ 
ten  as 

/,  -iknz  *  - iknz 

Ein  =  yE0e  ,  H\n  =  -xE0/r\0e  (1) 


The  electromagnetic  field  scattered  in  the  regions  Z  <  0  and  z  >  d  can  be  expressed  in  terms  of  the  Floquet  modes. 
Using  the  FEM  formulation,  the  weak  form  of  wave  equation  to  be  solved  for  the  electric  field  in  the  unit  cell  takes 
following  form  [5] 
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where  E  is  the  electric  field  in  the  unit  cell  volume  V  bounded  by  z  =  0  and  z  =  d  planes,  e  is  the  vector 

th  ^ 

modal  function  for  the  Floquet  modes,  Y p  is  the  modal  admittance  for  the  p  mode,  T  is  the  testing  function,  and 

,  5 2  are  the  surface  areas  at  the  z  —  0  and  z  —  d  planes,  respectively.  From  the  solution  of  equation  (2)  the 

reflection  and  transmission  coefficients  for  the  structure  shown  in  Figure  1(b)  can  be  determined. 

FDTD  Method:  It  is  often  required  to  estimate  the  reflection  and  transmission  of  EM  waves  through  meta-material 
over  multiple  frequencies.  The  FEM  approach  presented  above  becomes  computationally  intensive  when  required  to 
be  repeated  for  each  frequency  of  interest.  The  FDTD  method  provides  an  alternative  and  efficient  way  of  estimating 
reflection  and  transmission  through  meta-materials  over  multiple  frequencies.  To  perform  these  computations  using 
the  FDTD  approach  we  used  the  commercially  available  CST  Microwave  Studio  [6]  computational  tool.  To  validate 
the  FEM  code  and  the  numerical  results  obtained  using  the  CST  software,  a  dielectric  slab  embedded  with  four  split 
rings  and  capacitively  loaded  strip  dipole  as  shown  in  Figure  2(a)  is  considered.  The  reflection  and  transmission  coef¬ 
ficients  obtained  using  both  approaches  are  shown  in  Figure  2(b).  The  close  agreement  between  the  two  numerical 


results  confirms  validity  of  the  FEM  formulation  and  the  also  the  modeling  procedure  followed  while  using  the  CST 
Microwave  Studio. 

Using  these  techniques  the  reflection  and  transmission  characteristics  of  metallic  inclusions  of  various  shapes 
and  sizes  are  computed  and  results  will  be  shared  with  the  audience  at  the  time  of  presentation. 


Figure  2(a)  geometry  of  split  ring  and 
capacitive  loaded  strips. 
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Figure2(b)  Reflection,  transmission  characteristics  of 
split  ring/strip  geometry. 


Extraction  of  Electrical  Properties 

Following  the  Nicholson-Ross  [4]  method  to  estimate  electric  and  magnetic  properties  of  meta-materials,  it  is 
assumed  that  the  dielectric  slab  with  the  metallic  inclusions  is  replaced  by  a  homogeneous  material  slab  having  effec¬ 


tive  permittivity  |  e 
can  be  written  as 
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and  permeability  I  p  ^  J  .  The  reflection  and  transmission  coefficients  of  this  effective  slab 
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and  r  =  hi  /e  -1/1  hi  /e  +  1  .  The  problem  of  extraction  of  electric 
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properties  of  an  equivalent  slab  can  now  be  defined  as  determination  of  8f.  ,  \lr  using  (3)  from  the  knowledge  ol 


5ij,  S21,  d.  Defining  new  variables  V ^  =  S^+S0^,  V0  =  ^21_^11  anc*  ^  =  (1  _  3/9)/(V^  -  V9) ,  it 

can  be  shown  that 


r  =  Xljx2  -  1  and  t  =  (V1-r)/(l-V1r) 
The  steps  involved  in  extraction  process  are: 

(a)  determine  (r,  t)  using  (4), 
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(b)  express 
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(c)  with  proper  choice  of  n,  parameters  are  given  by  =  (( 1  +  r)/(  1  -  r))kr ,  e  =  (( 1  -  r)/(  1  +  r))kr 

Due  to  ambiguity  in  selecting  proper  value  of  n,  earlier  researchers  [7]  assuming  d  «  A.  had  simplified  the  calcula¬ 
tions  in  the  step  (b).  However,  the  assumption  d  «  X  limits  application  to  very  thin  slab.  In  this  work  we  have  devel- 


oped  a  procedure  to  track  correct  value  of  n  without  making  the  assumption.  The  proper  values  of  effective  electric 
properties  are  also  estimated  using  the  Genetic  Algorithm. 

Numerical  Results 

In  this  section  numerical  data  for  the  effective  parameters  of  a  slab  embedded  with  split  rings  arranged  in  a  peri¬ 
odic  fashion,  shown  in  Figure  3(a),  are  presented.  First,  the  unit  cell  including  a  single  layer  of  split  ring  is  modeled 
using  a  CAD  package.  The  transmission  and  reflection  coefficients  as  a  function  of  frequency  are  then  computed 
using  the  CST  Microwave  Studio  software.  Using  the  extraction  procedure  described  above  the  effective  parameters 
are  estimated  and  are  shown  in  Figure  3(b).  Using  the  extracted  parameters  the  reflection  and  transmission  coeffi¬ 
cients  are  calculated  and  compared  with  the  numerical  data  obtained  using  the  CST  Microwave  Studio. 


Thin  Metallic  Split  Rings 
Arranged  In  Periodic 
Fashion 


Composite  Slab  (Infinite  in  X-Y) 
Loaded  with  Metallic  Split 
Rings 


Figure  3(a)  Geometry  of  a  slab  loaded 
with  split  rings. 


Figure  3(b)  Effective  parameters  of  slab  loaded 
with  split  rings. 


Conclusions 

A  computational  electromagnetic  method  using  the  FEM  and  FDTD  approaches  have  been  presented  to  estimate 
the  reflection  and  transmission  through  a  meta-material  slab.  Nicholson-Ross  technique[4]  and  also  Genetic  Algo¬ 
rithm  procedure  in  conjunction  with  the  knowledge  of  reflection  and  transmission  coefficients  have  been  used  to 
extract  effective  permittivity  and  permeability  of  a  meta-material  slab.  The  reflection  and  transmission  coefficient  of 
actual  meta-material  slab  have  compared  well  with  the  numerical  data  obtained  using  the  extracted  effective  permit¬ 
tivity  and  permeability.  Due  to  page  limitation  the  numerical  results  for  a  metallic  inclusion  in  the  shape  of  split  ring 
only  has  been  shown  here.  The  metallic  inclusions  of  other  shapes  will  be  discussed  at  the  time  of  presentation. 
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